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Abstract 

We investigate percolation phenomena in multifractal objects that are built in a simple way. 
In these objects the multifractality comes directly from the geometrical tiling. We identify some 
differences between percolation in the proposed multifractals and in a regular lattice. There 
are basically two sources of these differences. The first is related with the coordination number 
that changes along the multifractal. The second comes from the way the weight of each cell in 
the multifractal affects the percolation cluster. We use many samples of finite size lattices and 
draw the histogram of percolating lattices against site occupation probability. Depending on a 
parameter characterizing the multifractal and the lattice size, the histogram can have two peaks. 
We observe that the percolation threshold for the multifractal is lower than one for the square 
lattice. We compute the fractal dimension of the percolating cluster and the critical exponent 
p. Despite the topological differences, we find that the percolation in a multifractal support is 
in the same universality class of the standard percolation. 

PACS numbers: 64.60.Ak 61.43.Hv 05.40.-a 91.60.-x 
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1- INTRODUCTION 



Percolation theory has been used in several fields such as chemistry, epidemics, 
science of materials, transport of fiuids in porous media, branched polymers, and econo- 
physics The original percolation model based on a square lattice has 

been extended to several kinds of regular and random lattices, to continuous media where 
the objects overlap in space, and to other complex systems BiEil In this work 
we generalize the percolation theory to cover an even broader range of complex systems. 
We devised an approach to investigate how percolation occurs in a support that is itself 
a multifractal. For this purpose we have constructed an easy assembling multifractal 
immersed in a 2D space. 

Our work is inspired in the modeling of geophysical natural objects that show multi- 
fractal properties 0, 0, JJ, 0|. The model can be applied to transport of fiuid in a 



multifractal porous media such as sedimentary strata. Oil reservoirs are possible candi- 
dates to be modeled in such a way since that the measurement of some physical quantities 
in well-logs show multifractal behavior 0, |^ . Despite the potential applications inspira- 
tion, this problem is important by itself in the scientific context. The study of percolation 
phenomena in multifractal lattices is relevant in Statistical Physics specially when the 
size of the blocks and their number of neighbors can vary. 

In order to make this analysis we create a multifractal object that can be used as a toy 
model and a laboratory for percolation theory. An important characteristic of this object 
is that its topologic properties (e.g. number of neighbors of each block) change over the 
object. In the reference jisl an algorithm that has some resemblance with ours is used. 
That multifractal is built from the partition of the square, but the object has a trivial 
topology. Besides, the object used in [3| is aleatory and ours is deterministic. Although 
both models present multifractality, our model has the following differences: it shows a 
non-trivial topology, we can determine analytically its spectrum of fractal dimensions, it 
generalizes the square lattice, and it shows simplicity in the construction. 

The multifractal object we have developed is a natural generalization of the regular 
square lattice once we consider the algorithmic point of view. The algorithm that gener- 
ates a square lattice with 2" x 2" cells starting with a square of fixed size is the following. 
We begin with a L x L square and we cut if in 4 identical pieces (cells). At each step 
all the cells are equally divided in 4 parts using vertical and horizontal segments. This 
process produces a lattice as a partition of the square. The multifractal we create is also 
a partition of the square, but the ratio we divide the cells is different from L/2. The 
parameter characterizing the multifractal, p, is related to the ratio of this division. 

What makes this problem new and appealing to physics is the following. The support 
of the percolation clusters is composed by subsets of different fractal dimensions. It is 
important to know how these different subsets are connected and how they participate 
in the conducting process. There are intriguing features in the network due to the fact 
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that all cells have rectangular shape but the area and the number of neighbors can vary, 
forming an exotic tilling. 

In section 2 we present the multifractal object that we construct to study percolation, 
and we analyze how its multifractal partition maps into the square lattice. In section 3 we 
expose the algorithm we use to estimate percolation threshold and derive the multifractal 
spectrum of the multifractal object. In section 4 we show numerical results and discuss 
the histograms of percolating lattices versus occupation probability. Finally in section 
5 we summarize the main differences between percolation in a regular lattice and in a 
multifractal support. 

2- THE MULTIFRACTAL OBJECT Qmf 

The central object of our analysis is a multifractal object that we call Qmf- Before 
defining it we enumerate some of its properties. 

• Qmf is a multifractal, it means, Qmf has an infinite number of fc-subsets each one 
with a distinct fractal dimension Dk- 

• It is possible to determine analytically the spectrum of all D^. 

• The sum of all the families of fc-subsets fills the square. This fact enables us to 
study its percolation properties using procedures similar to the ones applied in the 
site percolation in the square lattice. 

• The algorithm of construction of Qmf has just one parameter p. 

• For the special choice p = 1 the object Qmf degenerates into the square lattice. In 
this case we compare our results with the square lattice site percolation. 

• The object Qmf shows self-affinity or self-scaling depending on the region of the 
object. 

• Finally, the algorithm for construction of Qmf is simple and it is of easy implemen- 
tation in the computer. 

We define Qmf through the following algorithm. We start with a square of linear size 
L and choose a parameter < p < 1, where p = ^ for r and s integers. In the first step, 
n = 1, the square is cut in two pieces of area = and = by a vertical line 
(we use units). In other words, the square is cut according to a given p. This step is 
shown in figure [T] (a), in this figure we use p = ^ = |. 

In the second step, n = 2, we cut the two rectangles of figure [T] (a) by the same p, but 
using two horizontal lines as shown in figure ^ (b). This partition of the square generates 
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four rectangular blocks: the smallest one of area (j^)^, two of them of area (j^^p- and 
the largest one of area (j^)^, in the figure p > 0.5. 

The third step, n = 3, is shown in figure ^ (c) and the fourth step, n = 4, in (d). As 
observed in figure, at level n = 4 there are 2"^ blocks and the distribution of areas among 
the blocks follows the binomial law: 



1 + p 1 + P 1 + P 1 + P 1 + P ^l + p'^l + p' l + p 

We call the elements with the same area as a fc-set. In the case n = 4 we have five /c-sets. 

At step n the square has 2"~^ line segments, (n + 1) fc-sets and 2^ blocks. The partition 
of the area A = 1 (using units) of the square in different blocks follows the binomial 
rule: 

A = y C?{-^fi-^Y~^ = (1^)" = 1. (2) 
fc=0 1 + P 1+P 1 + P 

As n — >■ oo each fc-set (a subset made of cells of same area) determines a monofractal 

whose dimension we calculate in the next section. The ensemble of all fc-sets engenders 

the multifractal object Qmf- 

3- THE ALGORITHM OF PERCOLATION AND THE MULTIFRACTAL SPEC- 
TRUM 

In this section we show the algorithm used to study the percolation properties of 
Qmf and the analytic derivation of its spectrum of fractal dimensions. The estimation 
of the spectrum, D^, is performed using the box counting method 0] whose measure 
elements came from the percolation algorithm. 

The concept of the percolation algorithm for Qmf consists in mapping it into the 
square lattice. The square lattice should be large enough that each line segment of Qmf 
coincides with a line of the lattice. Therefore we consider that the square lattice is more 
finely divided than Qmf- In this way all blocks of the multifractal are composed by a 
finite number of cells of the square lattice. 

To explain the percolation algorithm we suppose that Qmf construction is at step n. 
We proceed the percolation algorithm by choosing at random one among the 2" blocks of 
Qmf- Once a block is chosen all the cells in the square lattice corresponding to this block 
are considered as occupied. Each time a block of Qmf is chosen the algorithm check if the 
occupied cells at the underlying lattice are connected in such a way to form an infinite 
percolation cluster. The algorithm to check the percolation is similar to the one used in 

Ei 
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For the estimation of the spectrum of an object X we use the box counting method 
j2^. The object X is immersed in the plane of real numbers 3?^, we use the trivial metric. 
Cover by just-touching square boxes of side length e. Let N{X) denote the number of 
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square cells of side length e which intersect X. If 

log NjX) log NjX) 

Dx = hm.^o — j 1 — = hrriL^^ — z — , (3) 

log ^ log L 

is finite, then Dx is the dimension of the X. 

In our case the object X is a fc-set. Remember that the fc-set corresponds to a set of 
rectangles of the same area. For a /c-set we have that A^^ is given by: 

where is the binomial coefficient that express the number of elements fc-type, and 
^{n-k) jg axea, of each element of this set. If the square is partitioned n times 
horizontal cuts and ^ vertical cuts) its size is L = (s + r)"/^. Combining all this 
information we have for the fractal dimension of each /c-set: 

_ log Cj s" r("-^) 

log [s + r)2 

In the r = s = 1 case all subsets of Qmf are composed by elements of the same area, 
square cells. In this way the object is formed by a single subset with dimension: 

D = lirrin^ J'^^ ^] = 2, (6) 

log (1 + 1)2 

This result is expected once in this particular case Qmj degenerates into the square lattice 
that has dimension 2. 

In figure 121 we show the picture of Qmf for p = |. We have used n = 12, in (a) the full 
object is shown, in (b) a zoom of an internal square of the object is illustrated. We have 
used the same color to indicate the elements of a same k-set. The funny tilling depicted 
in the figure is common to Qmf with different values of p. 

Figure El shows the spectrum of Dk for n = 400 calculated from equation (0). The 
use of increasing n does not change the shape of the curve, it only increases the number 
of k and makes the curve appear more dense. We use (s, r) = (2, 3) to illustrate the 
asymmetry of the distribution. The spectrum has a maximum close to p n. In this case 
|400 ~ 270. It means that the majority of mass of the multifractal is concentrated in the 
fc-sets around this value. The spectrum is typically asymmetric around its maximum. 
Only the case {s, r) = (1, 1) is symmetric and the asymmetry of Dk increases as — >■ 1, 
which is related to the area distribution among the blocks as we shall see in the next 
chapter. 



4- NUMERICAL SIMULATIONS 

In this chapter we focus our attention on the numerical results obtained from the 
algorithm exposed above. We are interested mainly in analyzing the percolating properties 
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of Qmf- Figure IH (a) shows the histogram of percolating lattices versus the occupation 
probability p. The area under the histogram is normalized to unity. We use n = 10 and 
average the results over 40000 samples. We consider that a lattice percolates when it 
percolates from top to bottom or from left to right. The histogram of percolating lattices 
at both directions is similar but slightly shifted to the right. This shift is common in 
percolation (see the reference [2^ for percolation in the square lattice). 

We show in figure 0] (a) the results of simulations for the following values of (s, r): 
(1, 1), which degenerates into the square lattice; and (2, 1), (4, 1), (6, 1) which correspond 
to truly multifractals. In this figure the histograms corresponding to (2, 1), (4, 1) and (6, 1) 
are shifted to the left compared to the histogram of (1, 1). The peak of the histogram for 
1, 1) corresponds, as expected, to the square lattice size percolation threshold, Pc = 0.597, 
3], since this case matches exactly the square lattice. The other values of pc are shown 
in Table /. 

The reason why Qmf, for diverse p, shows roughly the same Pc comes from the topology 
of the multifractal. The topology of a set of blocks is related with the coordination 
number, c, which is defined as the number of neighbors of each block Q]. Qmf has the 
property that c changes along the object and with p. However, we compute the average 
coordination number Cave- These results do not depend significantly on p, neither on n, the 
number of steps to build Qmf, which determines the number of blocks. The value found. 
Cave = 5.436, for the multifractal is close to the value of c of the triangular percolation 
problem which has c = 6 and whose analytic percolation threshold is pc = 0.5. The 
situation (s, r) = (1,1), the square lattice, shows trivially c = 4. Because the square 
lattice has a different c it configures a particular situation compared to other Qmf and it 
shows a different Pc as depicted in figure IH (a). 

In table / we show pc and the fractal dimension of the percolating cluster, df, for 
diverse p. We have done an average over 100000 samples and n = 16. The estimation of 
df is done by the relation M ~ L'^^ for the 'mass ' ,M , of the percolation cluster, that 
means, the area of the cluster measured in units of the underlying square lattice, and L, 
the size of the underlying lattice. Based on the values of rfj of table / we conclude that 
the percolation on a multifractal support (imbeded in two dimensions) belong to the same 
class of universality than the usual percolation in two- dimensions. The calculated value 
of df for the (6, 1) case is smaller compared to the others because of finite size effects. We 
discuss in detail this effect in the following paragraphs. 

Percolation shows critical phenomena and several scaling relations are observed. The 
critical exponent (3 is defined from the equation: 

i?i~(p,(L)-p,)^ (7) 

where Pc is the exact occupation probability value in contrast to Pc{L) which is the finite 
size value. The power-law (0) is verified for Pc{L) obtained from Rl. The numerical 
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estimation of (3 is based in equation (0) where Rl is a. key element of the analysis. For 
Qmf the probability Rl is not a well behaved function of p for low L as we shall see in 
the next paragraphs. Actually, Rl can show, depending on p, an inflection point at Pc in 
this regime. However, in the case where L ^ oo the scaling of {pc{L) — Pc) recovers an 
usual behavior. In this regime we find the same /5 characteristic of two dimension case, 
/3 = ^ = 0.13888. We checked in our simulations that for n = 18, (3 is around 5% of the 
exact value. The full set of values of (3 is in table /. 
Table I 



s,r) 


(1,1) 


(2,1) 


(3,1) 


(3,2) 


(4,1) 


(5,1) 


(6,1) 


Pc 


0.593 


0.527 


0.526 


0.526 


0.525 


0.525 


0.530 


df 


1.895 


1.900 


1.911 


1.890 


1.902 


1.929 


1.842 


P 


0.127 


0.128 


0.140 


0.141 


0.141 


0.118 


0.109 



It is worth to say that, despite small fluctuations in the values shown in the table, 
there is no tendency in the numbers. The conclusion we take from this data is that the 
errors are caused by finite size effects and low-number statistics. 

The dispersion of the histogram changes significatively with [s, r) as intuitively ex- 
pected. To illustrate the changing in the width of the histogram of a generic [s, r) multi- 
fractal we analyze the area of its blocks. At step n of the construction of Qmf the largest 
element has the area j^^^ and the smallest ^^^^ (using units). In this way the largest 
area ratio among blocks increases with (f)"'. As the occupation probability, entering in 
the percolation algorithm, is in general proportional to the area of the blocks we expect 
that the width of the histograms in figure 0] (a) increases with (^)". This increase in the 
dispersion is visualized clearly in the curves (2, 1) and (4, 1) of the figure. 

The most singular curve in figure 0] (a) is (6, 1) which shows clearly two peaks. We 
stress this point when we comment figure El Figure 111(b) uses the same data of figure (a), 
but instead of the histogram of percolating lattices we show the cumulative sum, Rl. As 
Rl is normalized, this parameter is also called the fraction of percolating lattices. As in 
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m (a) the case (s,r) = (1, 1), the square lattice, reproduces the results of literature 
In this situation the lattice size, L, is L = {s + r)^° = 1024. For this special case the 
number of blocks is equal to the number of unit boxes covering the surface. The double 
peak case (s,r) = (6, 1) shows an infiection point in the graphic of Rl versus p. In the 
following figure we explore in detail this point. 

The most noticeable signature of percolation in the multifractal Qmf is the double peak 
observed for (s, r) = (6, 1) in figure El In this figure the histograms of percolating lattices 
versus p is plotted for diverse n as indicated in the figure. The distance between the peaks 
decreases as n increases. This picture indicates that the double peak is a phenomenon 
that is relevant for percolation in the multifractal, when p is low, in the finite lattice size 
condition used in the simulation. From an analytic point of view the curve (6, 1) in figure 
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Elis different from curve (1, 1). In curve (6, 1) there are three extrema points while in the 
(1, 1) case the curve shows a single maximum point. We conjecture that in the limit of 
n — * oo these three points coalesce into a single one and all the curves show a similar 
behavior. 

The two peaks in the histogram come from the huge difference among the area of 
the blocks of Qmf- For large (^)"' the area difference is so accentuated that we model 
the histogram of percolating lattices as a bimodal statistics. In the case of the largest 
block be chosen the multifractal easily percolates compared with the opposite possibility. 
To estimate the effect of the largest area block on the statistic we make Table //. The 
difference between the first peak at pi and the second one at p2 is Apmax- In Table // we 
compare Apmax with the fraction of the largest block over the total square area (^)"- 
This comparison is made for different steps in the construction of the multifractal n, as 
n increases the area difference decreases as well as the distance between peaks. Table // 
shows a good agreement between both values, we conclude that the bimodal statistic is 
caused by the huge mass of the largest block. 

Table II 

n 8 10 12 14 16 18 

Ap^ax 0.29 0.22 0.15 0.11 0.070 0.040 

0.291 0.211 0.157 0.115 0.084 0.062 

We notice, however, that the agreement between Apmax and (7^)"' decreases as n 
increases. We interpret the disagreement between the bimodal statistics hypothesis and 
the numerics for high n as the limit of the hypothesis. Actually, the largest block is not 
the only one that produces anisotropy in the multifractal, and as n increases this fact 
becomes more accentuated. For small n the large block can be taken as the main factor of 
anisotropy, and the bimodal statistics apply. Large n implies, however, truly multifractals 
and a more complex statistics should be used to treat the problem. 

5- CONCLUSION 

In this work we develop a multifractal object, Qmf to study percolation in a mul- 
tifractal support. Besides of Qmf being a multifractal, it shows several interesting prop- 
erties. The sum of all its fractal subsets fills a square and it is possible to determine the 
spectrum of its fractal dimensions. In addition, the algorithm that generates Qmf has 
only one free parameter p, and in the p = 1 case Qmf becomes the square lattice. 

We observe that percolation in a multifractal presents different features from percola- 
tion in a regular lattice. There are two reasons for that: the heterogeneous distribution 
of weight (area) among the blocks and the variation of the coordination number of the 
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topologic structure. The weight of each block in a multifractal counts diversely in the 
mass of the infinite percolating cluster. The difference in weight of the blocks changes the 
dispersion of the histogram of percolating lattices. The phenomenon of two peaks appear- 
ing in the histogram is also connected with the weight difference. We model the distance 
between peaks using a bimodal statistics. In the limit of n — >^ oo all the histograms of 
multifractals seems to collapse into a single curve. 

For all cases in which p 7^ 1 the multifractal Qmf shows a coordination number (number 
of neighbors of each block) that changes along the object. The average coordination 
number of Qmf is around 5.436. In contrast, the situation (s,r) = (1,1) (the square 
lattice particular case) has a coordination number constant and equal to 4. This suggests 
that the case p 7^ 1 represents a break of the symmetry of the system. In this sense the 
coordination number (topology) is much more complex for Qmf than for a regular lattice. 
Despite these differences, we have done numerical estimations of the fractal dimension of 
the percolating cluster in the multifractal, obtaining values which are around 1.89, the 
same dimension found for the incipient percolation cluster in a two-dimensional regular 
lattice. The numerical simulation of the (3 critical exponent, also show the same value 
of two-dimensional regular case, and point to the same conclusion that we have regular 
percolation. 
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FIGURE LEGENDS 

FIG. 1: The four initial steps in the formation of Qmf- In (a) the vertical line cutting the 
square in two pieces of area ratio p. Two horizontal lines sharing the rectangles by the same 
ratio are depicted in (b). The third step is indicated in (c) and the fourth in (d). At each step 
the areas of the respective blocks are shown in figure. 



FIG. 2: The figure shows the three views of the multifractal Qmj, it is used n = 12, (s, r) = (3, 2). 
In (a) we have the original picture, figure (b) is a zoom of the square indicated in (a). 



FIG. 3: The spectrum of fractal dimensions Dk of Qmf for n = 400 and (p, q) = (3, 2). 
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FIG. 4: In (a) is depicted the histogram of percolation lattices versus the occupation probability 
p for the cases (s,r) = (1, 1), (2, 1), (4, 1), and (6, 1). The areas under the curve are normalized 
to unity. For the same (s,r) is shown in (b) the graphic of the fraction of percolation lattices Rl 
versus p. It is used 40000 lattices to make the average. 



FIG. 5: The histogram of percolated lattices versus the occupation probability p for several 
values lattice size. The graphic shows the double peak approaching each other as n increases, 
in the figure (s,r) = (1,6) and 8 < n < 18. It is used 40000 lattices to make the average. 
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